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We  report  experimental  observations  of  synchronization  among  combustion  variations  in  different  cyl¬ 
inders  at  fuel-lean  conditions  in  an  eight-cylinder  spark  ignition  engine.  Our  results  appear  to  confirm  that 
synchronization  readily  occurs  and  that  it  becomes  stronger  as  the  overall  equivalence  ratio  is  reduced 
from  stoichiometric.  It  also  appears  that  the  onset  of  synchronization  is  associated  with  bifurcation  insta¬ 
bilities  reported  previously  for  combustion  in  single  cylinders.  We  use  both  cross-correlation  and  symbolic 
time  series  analysis  to  quantify  the  apparent  relationships  between  pairs  of  cylinders  and  multicylinder 
groups.  Extension  of  a  simple  dynamic  model  for  single-cylinder  combustion  variations  to  the  multicylinder 
case  appears  to  agree  with  the  observations  and  provides  a  basis  for  further  studies.  The  occurrence  of 
significant  cylinder-to-cylinder  synchronization  may  have  significant  implications  for  engine  diagnostics  and 
control. 


Introduction 

Fuel-lean  operation  of  spark  ignition  engines  is  of 
growing  importance  for  reducing  emissions  and  im¬ 
proving  fuel  efficiency.  One  of  the  major  practical 
constraints  to  lean  operation  has  been  the  onset  of 
combustion  instabilities  near  the  lean  limit,  resulting 
in  an  unacceptable  level  of  misfires  and  partial 
burns.  Thus,  there  is  considerable  interest  in  under¬ 
standing  the  fundamental  causes  of  lean  combustion 
instability  and  how  the  negative  effects  can  be  min¬ 
imized. 

Recent  studies  of  cycle-to-cycle  combustion  vari¬ 
ations  within  a  single  cylinder  have  shown  that  the 
behavior  under  lean  fueling  conditions  is  consistent 
with  the  onset  of  a  noisy  period-doubling  bifurcation 
sequence  [1-5].  Such  instabilities  are  a  common  fea¬ 
ture  of  nonlinear  dynamical  systems  in  which  future 
behavior  is  highly  sensitive  to  small  changes  in  past 
conditions.  In  the  case  of  lean  premixed  engine  com¬ 
bustion,  residual  fuel  and  air  have  a  highly  nonlinear 
effect  on  subsequent  ignition  and  flame  propagation 
in  succeeding  combustion  strokes.  As  has  been  pre¬ 
viously  demonstrated,  the  resulting  cycle-to-cycle 
dynamics  can  be  simulated  using  a  simple  mapping 
function  that  predicts  the  outcome  of  future  com¬ 
bustion  events  (within  a  given  cylinder)  based  on  en¬ 
gine  operating  parameters  and  the  outcome  of  past 


combustion  events  (in  the  same  cylinder).  The  ef¬ 
fects  of  higher-order  phenomena  such  as  in-cylinder 
fluid  turbulence  and  vaporization  of  fuel  droplets  are 
represented  as  stochastic  perturbations  of  engine  pa¬ 
rameters  such  as  residual  fraction  and  as-injected 
equivalence  ratio.  The  purpose  of  such  a  hybrid  low- 
order  model  is  to  capture  the  essential  statistical  fea¬ 
tures  of  the  combustion  instability  without  being 
overwhelmed  by  details.  Thus,  one  sacrifices  infor¬ 
mation  about  mixing  and  fluid  turbulence  for  rapid 
prediction  and  the  evaluation  of  large  numbers  (e.g., 
thousands)  of  cycles.  This  type  of  model  is  expected 
to  be  most  useful  for  understanding  general  trends 
(e.g.,  for  engine  diagnostics)  or  for  adaptive  feedback 
control.  Practical  control  applications  can  be  seen  in 
Refs.  [6,7],  which  describe  how  cycle-resolved  feed¬ 
back  can  be  used  to  counter  combustion  instabilities. 

In  this  study,  we  were  interested  in  determining  if 
lean  combustion  instability  is  also  associated  with 
cylinder-to-cylinder  coupling.  Specifically,  such  cou¬ 
pling  can  lead  to  a  condition  known  as  synchroni¬ 
zation,  in  which  combustion  events  in  one  or  more 
cylinders  develop  some  type  of  predictable  pattern 
relative  to  each  other.  This  investigation  differs  from 
earlier  studies  [I-5,8,9]  in  that  we  are  investigating 
how  cylinders  behave  collectively  rather  than  indi¬ 
vidually.  If  synchronization  is  confirmed,  our  ulti¬ 
mate  objective  is  to  develop  real-time  diagnostics 
and  controls  for  this  collective  behavior. 
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Fig.  1.  Cylinder  and  intake  manifold  configuration.  Cyl¬ 
inder  location  is  denoted  by  C  for  physical  location  and  E 
for  location  in  firing  order. 

Synchronization  is  a  general  property  of  multi¬ 
element  nonlinear  systems,  and  many  examples  have 
been  reported  [10-14].  Perhaps  the  most  famous  ex¬ 
ample  is  the  case  of  the  synchronizing  clocks  re¬ 
ported  by  Christian  Huygens  in  1673  [13].  Huygens 
demonstrated  that  clocks  mounted  on  a  common 
support  were  synchronized  through  small  mechani¬ 
cal  vibrations  that  were  exchanged  between  them. 
Similar  behavior  has  been  reported  in  cellular  flames 
[14]  and  cardiorespiratory  systems  [12]. 

Our  approach  included  experimental  combustion 
measurements  from  a  multicylinder  engine,  analysis 
of  the  resulting  data,  and  comparison  of  our  exper¬ 
imental  observations  with  predictions  from  a  low- 
order  model.  In  the  following  sections,  we  describe 
our  experimental  setup,  data  analysis  procedures, 
and  the  proposed  model. 


Experimental  Apparatus  and  Procedures 

The  experiments  discussed  here  were  conducted 
on  a  1999  4.6  L,  two-valve  V-8  Ford  Expedition  en¬ 
gine.  The  nominal  operating  conditions  for  the  ex¬ 
periments  were  1200  rpm,  25  N-m  brake  torque, 
and  20°  BTDC  (before  top  dead  center)  spark  tim¬ 
ing.  The  engine  had  production  port  fuel  injection 
with  a  bore  of  90.2  mm,  a  stroke  of  90.0  mm,  a  com¬ 
pression  ratio  of  9.0,  and  a  valve  overlap  of  27  de¬ 
grees.  Engine  speed  was  controlled  with  an  absorb¬ 
ing/motoring  General  Electric  (GE)  dynamometer. 

In-cylinder  pressure  measurements  were  re¬ 
corded  once  per  crank  angle  degree  at  steady-state 
conditions  as  equivalence  ratio  was  varied  over  five 
equivalence  ratios  from  stoichiometric  to  very  lean 
fueling  conditions.  The  firing  order  for  the  cylinders 


was  1-3-7-2-6-5-4-8  (see  Eig.  1  for  the  cylinder  lay¬ 
out).  Eeedback  control  was  engaged  to  achieve  an 
operating  condition;  once  the  condition  was 
achieved,  the  feedback  controllers  were  shut  off,  and 
the  engine  was  run  open-loop.  The  only  controller 
operating  during  data  acquisition  was  the  dynamom¬ 
eter  speed  controller,  which  actuated  dynamometer 
torque  to  maintain  constant  engine  speed.  This  strat¬ 
egy  ensured  that  the  engine  combustion  dynamics 
were  minimally  influenced  by  controllers. 

The  memory  buffers  of  the  data  acquisition  system 
limited  contiguous  data  sets  to  approximately  354 
engine  cycles  (for  all  eight  cylinders).  Repeat  sets  of 
contiguous  data  blocks  were  collected  after  a  2  or  3 
min  pause.  For  the  results  described  below,  three 
contiguous  data  blocks  were  acquired  at  equivalence 
ratios  of  1.00,  0.83,  0.71,  0.66,  and  0.59.  All  experi¬ 
ments  were  then  repeated  a  second  time  to  evaluate 
reproducibility.  Net  heat  release  values  for  each  cy¬ 
cle  were  calculated  from  the  cylinder  pressure  data 
utilizing  a  variant  of  the  Rassweiler  and  Withrow  in¬ 
tegration  method  [15]. 


Experimental  Results 

Figure  2a  and  b  illustrates  typical  heat  release  pat¬ 
terns  in  all  eight  cylinders  at  stoichiometric  and  very 
lean  fueling  conditions.  These  plots  depict  pairs  of 
sequential  heat  release  values  for  each  cylinder  and 
illustrate  how  each  combustion  event  relates  to  its 
successor  in  time.  We  note  that  at  near  stoichiomet¬ 
ric  fueling  (Fig.  2a),  combustion  variations  are  very 
small,  producing  a  focused  cluster  of  points  around 
the  average  value.  Detailed  analysis  of  these  small 
variations  shows  that  they  have  a  Gaussian  distribu¬ 
tion  and  appear  to  occur  randomly  in  time. 

At  very  lean  fueling  (Fig.  2b),  the  picture  changes 
considerably.  As  shown  in  Fig.  2,  we  observe  a  pat¬ 
tern  of  variations  that  is  both  far  from  Gaussian  dis¬ 
tributed  and  clearly  not  random.  The  detailed  pat¬ 
terns  vary  from  cylinder  to  cylinder,  but  all  cylinders 
exhibit  the  characteristic  hooked  shape  reported 
previously  for  single-cylinder  studies  [3,5] .  Our  anal¬ 
yses  for  each  individual  cylinder  also  confirm  that  the 
statistical  properties  of  the  measurements  are  con¬ 
sistent  with  the  noisy  bifurcation  model  of  Daw  et 
al.  [!].  However,  for  this  multicylinder  engine  there 
is  one  complication;  each  cylinder  appears  to  pro¬ 
gress  through  the  bifurcation  sequence  at  a  different 
rate  as  fueling  is  leaned.  Thus,  at  each  fueling  con¬ 
dition,  each  cylinder  is  bifurcated  to  a  greater  or 
lesser  extent  and  exhibits  a  distinctive  pattern  of 
combustion  variations.  We  interpret  these  cylinder- 
to-cylinder  differences  as  resulting  from  non-uni¬ 
formities  in  the  distribution  of  fuel  and  air  (espe¬ 
cially).  Such  non-uniformities  are  well  known  to 
occur  and  pose  major  challenges  to  engine  designers 
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Fig.  2.  Heat  release  return  maps  for  all  eight  cylinders  at  (a)  stoichiometric  and  (b)  very  lean  fueling  conditions,  </>  = 
0.66.  Numbers  in  (a)  refer  to  cylinder  location. 


Simultaneous  measurements  of  all  eight  cylinders 
allowed  us  to  investigate  the  possibility  that  com¬ 
bustion  variations  in  different  cylinders  might  be  re¬ 
lated.  We  began  by  looking  for  cross-correlations  be¬ 
tween  cylinders,  as  depicted  in  Fig.  3.  In  Fig.  3,  we 
illustrate  example  behavior  over  1000  engine  cycles 
for  selected  cylinder  pairs  as  fueling  was  leaned. 
Note  that  for  one  pair  (F-0,  F-5),  a  significant  posi¬ 
tive  correlation  developed  when  the  as-injected 
equivalence  ratio  was  adjusted  to  0.66.  Likewise,  a 
significant  anticorrelation  appeared  simultaneously 
for  the  cylinder  pair  (F-0,  F-1),  while  there  appeared 
to  be  no  significant  correlation  for  the  pair  (F-0,  F- 
4).  Although  all  possible  pairs  are  not  shown  in  Fig. 
3,  we  observed  many  cylinder  pairs  with  cross-cor¬ 
relations  exceeding  0.2  under  lean  conditions.  We 
can  summarize  our  observations  from  such  cross¬ 
correlations  as  follows: 

•  At  stoichiometric  fueling,  we  did  not  observe  sta¬ 
tistically  significant  cross-correlations  between  any 


cylinder  pairs,  implying  that  combustion  variations 
in  each  cylinder  were  effectively  independent. 

•  As  fueling  was  leaned,  significant  cross-correla¬ 
tions  (both  positive  and  negative)  developed  be¬ 
tween  combustion  variations  in  some  cylinder 
pairs. 

•  The  distribution  of  positive  and  negative  cross-cor¬ 
relations  among  the  different  cylinder  pairs  was 
complex  and  did  not  have  any  obvious  relationship 
to  cylinder  location  or  firing  order.  We  observed 
strong  correlations  between  pairs  fed  by  the  same 
or  opposite  sides  of  the  intake  plenum,  pairs  on 
the  same  or  opposing  sides  of  the  engine,  and  pairs 
close  or  far  apart  in  firing  order.  Note  that  the  pairs 
(F-0,  F-5)  and  (F-3,  F-4)  have  similar  physical  cyl¬ 
inder  locations  with  respect  to  each  other,  but  the 
former  pair  shows  strong  correlation  under  lean 
conditions  while  the  latter  pair  shows  no  significant 
correlation. 

•  The  distribution  of  lean-fueling  cross-correlations 
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Fig.  3.  Cross-correlation  beUveen  selected  cylinder  pairs 
as  a  function  of  equivalence  ratio.  Significant  positive  and 
negative  cross-correlations  occurred  at  an  equivalence  ratio 
of  0.66. 

shifted  over  time.  Strong  correlations  (or  anticor¬ 
relations)  developed  between  specific  pairs,  per¬ 
sisted  for  a  significant  interval  (e.g.,  tens  to  hun¬ 
dreds  of  cycles),  and  then  abruptly  reduced  in 
magnitude. 

•  At  extremely  lean  fueling,  combustion  became 
very  poor  in  all  of  the  cylinders,  and  the  cross¬ 
correlations  between  cylinders  diminished. 

The  lack  of  significant  cross-correlations  at  stoi¬ 
chiometric  fueling  implies  that  any  potential  cou¬ 
pling  between  cylinders  was  initially  absent  or  at 
least  had  no  impact  on  combustion.  Likewise,  the 
occurrence  of  apparently  significant  cylinder-to-cyl- 
inder  correlations  at  lean  fueling  implies  that  some 
sort  of  coupling  mechanism  was  activated,  allowing 
some  of  the  cylinders  to  synchronize  with  each  other. 
One  possible  explanation  is  that  uncharacteristically 
large  or  small  pressure  oscillations  are  induced  in  the 
exhaust  and  intake  manifolds  when  the  combustion 
oscillations  in  individual  cylinders  become  promi¬ 
nent  (i.e.,  when  partial  or  total  misfires  occur).  Such 
changes  in  the  pressure  oscillations  could  impact  the 
subsequent  intake  and  exhaust  strokes  of  other  cyl¬ 
inders,  thereby  imposing  additional  variations  in  the 
effective  in-cylinder  equivalence  ratio  and  stimulat¬ 
ing  the  combustion  to  change.  The  apparent  varia¬ 
tions  in  bifurcation  level  from  cylinder  to  cylinder 
would  be  expected  to  complicate  the  coupling  inter¬ 
actions  because  the  stronger  oscillating  cylinders 
would  tend  to  drive  those  with  weaker  oscillations. 
While  we  did  not  measure  the  manifold  pressure 
with  sufficient  accuracy  to  confirm  this  scenario  in 
the  current  experiments,  the  impact  of  manifold 
pressure  waves  (acoustics)  on  engine  combustion  is 
a  widely  accepted  phenomenon  [15]. 


To  further  resolve  the  dynamic  relationships  be¬ 
tween  cylinders,  it  is  useful  to  apply  more  sophisti¬ 
cated  analyses.  One  method  we  find  to  be  particu¬ 
larly  useful  is  symbol  sequence  analysis  [2,3].  Briefly, 
symbolization  transforms  an  original  time  series  of 
measurements  into  a  sequence  of  discretized  sym¬ 
bols.  Such  a  transformation  makes  it  possible  to 
readily  detect  and  characterize  the  occurrence  of 
non-random  patterns,  even  when  there  is  consider¬ 
able  experimental  noise. 

In  the  case  of  multicylinder  combustion  measure¬ 
ments,  we  are  interested  in  observing  non-random 
patterns  over  both  space  and  time.  One  useful  sym¬ 
bolization  transforms  the  observed  heat  release  val¬ 
ues  into  a  binary  data  stream  (0  or  1)  depending  on 
where  they  fall  relative  to  the  median  value.  For 
each  cylinder,  we  use  the  median  value  for  that  cyl¬ 
inder  so  that  biases  in  the  nominal  fueling  levels  are 
discounted.  Within  the  symbolized  data,  we  evaluate 
interactions  between  specific  cylinder  pairs  by  con¬ 
structing  symbol  words  (i.e.,  short  symbol  se¬ 
quences)  according  to  the  following  definition: 

S2(j)  =  [x(j  -  1),  yij  -  1),  xij),  y{j), 

x{j  +  1),  yij  +  1)]  (1) 

where  S2(j)  is  the  symbolic  state  of  the  cylinder  pair 
at  engine  cycle j,  and  xij)  and  yij)  are  the  symbolic 
heat  release  values  at  cyclej  for  the  first  and  second 
cylinder  pair  members,  respectively.  For  conven¬ 
ience,  we  refer  to  each  symbol  sequence  according 
to  its  index  number  defined  by: 

Index(S2)  =  32%(j  —  1)  -I-  16ij(j  —  1) 

-I  8xij)  +  4yij) 

+  2x(j  -6  1)  -1  yij  +  1)  (2) 

This  scheme  allows  us  to  uniquely  refer  to  any  par¬ 
ticular  pattern  with  a  single  number  between  0  and 
63. 

Figure  4  illustrates  how  symbol  sequence  histo¬ 
grams  can  be  used  to  depict  the  correlated  and  un¬ 
correlated  cylinder  pairs  shown  previously  in  Fig.  3. 
On  the  abscissa  we  plot  the  sequence  index  number, 
and  on  the  ordinate  we  plot  the  relative  frequency 
at  which  each  possible  sequence  is  observed.  At  stoi¬ 
chiometric  fueling,  all  sequences  are  equally  prob¬ 
able,  and  the  observed  frequencies  lie  close  to  the 
heavy  dashed  line  in  Fig.  4.  When  combustion  os¬ 
cillations  develop  at  lean  fueling,  peaks  become 
prominent  at  sequence  index  values  12,  25,  38,  and 
51,  corresponding  respectively  to  the  sequences 
[0,0,1, 1,0,0],  [0,1,1,0,0,1],  [1,0,0,1,1,0],  and 

[1, 1,0,0, 1,1].  For  uncorrelated  behavior  (F-0,  F-4), 
the  heights  of  these  target  peaks  are  all  similar;  that 
is,  there  is  no  favored  relationship  between  the  cyl¬ 
inder  phases.  Thus,  even  when  there  are  consider¬ 
able  oscillations  in  each  cylinder,  the  lack  of  coupling 
means  that  we  are  just  as  likely  to  observe  in-phase 
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Sequence  Index 

Fig.  4.  Example  symbol  sequence  histograms  for  se¬ 
lected  cylinder  pairs  at  very  lean  fueling  conditions,  (fi  = 
0.66.  Solid  symbols  and  solid  line  correspond  to  cylinder 
pair  (F-0,  F-5)  which  had  a  significant  positive  correlation. 
Open  symbols  and  broken  line  correspond  to  pair  (F-O,  F- 
4)  which  had  no  significant  correlation.  The  heavy  dashed 
line  represents  the  case  where  all  sequences  are  equally 
probable,  which  was  approached  at  stoichiometric  fueling. 

or  out-of-phase  behavior.  When  there  is  significant 
correlation  (F-0,  F-5),  the  relative  heights  of  peaks 
12  and  51  are  higher  than  peaks  25  and  38,  indicating 
that  the  combustion  variations  in  both  cylinders  tend 
to  remain  in  phase  more  frequently  than  out  of 
phase.  Peaks  25  and  38  are  favored  for  the  same 
reason  when  the  cylinders  become  anticorrelated.  In 
all  cases,  the  mean  height  of  the  four  peaks  above 
background  levels  indicates  the  degree  of  joint  bi¬ 
furcation  appearing  in  both  cylinders.  It  is  thus  pos¬ 
sible  to  use  symbol  sequence  histograms  to  simul¬ 
taneously  evaluate  both  the  degree  of  bifurcation 
and  possible  coupling. 

Figure  5a  and  b  illustrates  how  symbol  synchro¬ 
grams  can  be  constructed  to  observe  bifurcation  and 
synchronization  over  time.  We  plot  the  observed  se¬ 
quence  index  against  engine  cycle  number,  for  ex¬ 
ample,  uncorrelated  and  anticorrelated  cylinder 
pairs,  respectively.  For  the  uncorrelated  pair,  we  ob¬ 
serve  that  all  four  target  peaks  occur  frequently,  but 
these  occurrences  are  sporadic  and  distributed 
evenly  over  time.  Thus,  both  cylinders  are  oscillating 
significantly,  but  the  oscillations  are  not  related  to 
each  other  in  any  consistent  way.  Conversely,  for  the 
anticorrelated  pair,  we  observe  significantly  long  un¬ 
interrupted  periods  (episodes)  of  relatively  strong 
anticorrelated  behavior  (peaks  25  and  38).  Such  ep¬ 
isodic  correlation  is  suggestive  of  noisy  synchroni¬ 
zation. 

Synchronization  among  groups  of  cylinders  larger 
than  two  can  also  be  studied  using  a  modified  symbol 


Cycle  Number 
(a) 


Cycle  Number 
(b) 


Fig.  5.  Symbol  synchrograms  for  (a)  uncorrelated  and 
(b)  anticorrelated  cylinder  pairs  at  very  lean  fueling  con¬ 
ditions,  </)  =  0.66.  Dark  bands  of  repeated  symbol  se¬ 
quences  reveal  combustion  oscillations.  A  persistent  bias  in 
these  bands  indicates  apparent  synchronization. 

sequence  definition.  Specifically,  we  define  a  three- 
cylinder  word  as 

S3(j)  =  [x(j  -  1),  y{j  -  1),  z(j  -  1), 

*(j).  1/(7).  2(j)]  (3) 

where  x{  j),  y{  j),  and  z(j)  are  symbolic  heat  release 
values  for  cylinders  x,  y,  and  z  measured  in  engine 
cyclej.  The  corresponding  sequence  index  is 

Index(S3)  =  32x(j  —  1)  -k  16y(j  —  1) 

+  8z{j  —  1)  -f  4x(j) 

+  %(;')  +  2(7)  (4) 

Fig.  6  is  an  example  synchrogram  based  on  the  above 
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Fig.  6.  Symbol  synchrogram  for  three  cylinders  (F-0,  F- 
4,  F-6)  at  very  lean  fueling  conditions,  </>  =  0.66.  We  ob¬ 
serve  a  persistent  bias  in  the  relationship  of  the  three  cyl¬ 
inders  as  a  group. 


Fig.  7.  Symbol  synchrogram  for  anticorrelated  model 
simulation.  Persistent  bands  reveal  the  anticorrelated  cou¬ 
pling  deliberately  added  to  the  model.  The  observed  pat¬ 
terns  are  qualitatively  similar  to  the  experimental  obser¬ 
vations  for  anticorrelated  cylinder  pairs. 

sequence  definition  applied  to  two  cylinder  pairs 
with  a  common  member.  In  this  case,  there  appears 
to  be  an  episodic  interaction  among  all  three  cylin¬ 
ders.  Patterns  similar  to  Fig.  6  were  also  observed 
for  other  three-cylinder  groups  in  our  engine. 

Proposed  Model  for  Multicylinder 
Synchronization 

Based  on  the  above  observations,  we  conjecture 
that  a  model  for  simulating  the  observed  behavior 


can  be  constructed  by  adapting  the  Daw  et  al.  [1] 
model.  Specifically,  we  suggest  that  the  random 
noise  term  used  for  representing  unmodeled  engine 
perturbations  should  be  modified  to  consist  of  a  ran¬ 
dom  component  and  a  non-random  component  that 
reflects  intercylinder  coupling.  From  a  dynamical 
systems  perspective,  this  implies  that  multicylinder 
engines  are  more  accurately  described  as  noisy  cou¬ 
pled  map  lattices,  where  each  cylinder  is  a  nonlinear 
mapping  function  that  is  coupled  in  varying  degrees 
to  the  others.  Such  a  model  might  be  useful  for  es¬ 
timating  the  average  coupling  between  cylinders 
(e.g.,  for  data  fitting).  It  is  also  conceivable  that  adap¬ 
tive  versions  of  this  model  (e.g.,  implemented  as  a 
neural  network)  might  be  used  for  on-board  diag¬ 
nostics  and  control. 

To  construct  such  a  model,  we  iterated  two  repli¬ 
cates  of  the  Daw  et  al.  map,  representing  two  dif¬ 
ferent  cylinders  in  a  noisy  bifurcated  state.  We  fur¬ 
ther  assumed  that  the  dominant  noisy  parametric 
inputs  to  both  cylinders  were  variations  in  the  as- 
injected  equivalence  ratio.  These  equivalence  ratio 
variations  were  assumed  to  be  the  sum  of  Gaussian 
variations  (independent  for  each  cylinder)  and  a 
small  oscillating  component  that  was  either  corre¬ 
lated  or  anticorrelated  between  the  cylinders.  The 
oscillating  component  was  intended  to  represent 
manifold  pressure  waves  that  might  effectively  cou¬ 
ple  the  cylinders  [16].  By  varying  the  relative  con¬ 
tribution  of  the  oscillatory  component,  we  studied 
the  predicted  effect  of  coupling.  The  level  of  total 
noise  input  to  each  cylinder  was  kept  constant. 

Example  results  of  the  above  model  are  shown  in 
synchrogram  form  in  Fig.  7.  An  anticorrelated  oscil¬ 
lating  noise  component  was  input  at  5%  (i.e.,  the 
percentage  of  noise  variance  contributed  by  the  os¬ 
cillating  component).  At  these  conditions,  we  ob¬ 
served  apparent  synchronization  behavior  that  was 
very  similar  to  some  of  our  engine  measurements  for 
anticorrelated  cylinder  pairs  (see  Fig.  5b).  Other 
similarities  observed  in  the  model  include: 

•For  both  the  model  and  experiment,  the  onset  of 
significant  cross-correlations  coincides  with  bifur¬ 
cations. 

•  Both  the  model  and  experiment  exhibit  transient 
correlation  periods  of  similar  length  (e.g.,  50-200 
cycles). 

•  The  statistical  frequencies  of  the  correlation  epi¬ 
sodes  are  similar  (e.g.,  interacting  cylinders  can  be 
synchronized  more  than  20%  of  the  time  once  bi¬ 
furcation  occurs). 

•  The  model  predicts  episodic  correlation  of  multi¬ 
ple  cylinder  groups  (more  than  two  at  a  time), 
which  we  have  experimentally  observed  (see  Fig. 
6). 

Thus,  it  appears  that  the  model  predicts  behavior 
that  is  at  least  qualitatively  correct. 
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Summary  and  Conclusions 

Our  experimental  observations  indicate  significant 
cylinder-to-cylinder  synchronization  where  the  onset 
of  synchronization  corresponds  to  the  development 
of  bifurcations  in  individual  cylinders.  Similar  be¬ 
havior  has  been  recently  confirmed  on  a  significantly 
different  engine  with  a  different  intake  manifold 
[17].  A  simple  mapping  model  with  a  shared  noise 
component  to  reflect  intercylinder  coupling  exhibits 
many  of  the  characteristic  features  observed  exper¬ 
imentally.  More  experimental  data  are  needed  to  re¬ 
fine  the  picture  of  the  possible  cylinder-to-cylinder 
interaction  patterns.  Future  studies  should  address 
the  physical  mechanisms  responsible  for  intercylin¬ 
der  coupling.  The  observation  of  the  reported  pat¬ 
terns  here  in  two  different  multicylinder  spark  ig¬ 
nition  engines  could  have  important  implications  for 
engine  design  and  control. 
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